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ABSTRACT 

Long ago, McVittie had found a class of solutions which can be 
thought of as Schwarzschild black holes in an FRW universe. In 
recent years they have been studied extensively and generalised 
to charged and uncharged black holes in D > 4 dimensions also. 
Here, assuming an ansatz similar to McVittie's, we present solu- 
tions for uncharged branes which can be thought of as branes in 
a time dependent universe. We consider their application to the 
brane antibrane decay process, also referred to as tachyon con- 
densation, and discuss the necessary generalisations required for 
our ansatz to describe such a process. 
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1. Static brane solutions in string/M theory are well known [1]. It is 
of obvious interest to study time dependent solutions also. In recent years 
there have been numerous studies of strings and branes in a variety of time 
dependent backgrounds [2]. 

Here we consider a different class of solutions, of the type McVittie found 
long ago for four dimensional Schwarzschild black holes [3]. These McVittie 
solutions are inhomogeneous and time dependent, and are thought of as 
describing black holes in an FRW universe with perfect fluid matter. These 
solutions and their properties have been discussed extensively in [4, 6, 5] and 
are also generalised to both charged and uncharged black holes in D > 4 
dimensions [7, 8, 9]. 

In this paper, we study brane solutions in a time dependent universe by 
generalising the McVittie ansatz suitably. We consider only the uncharged 
branes since the appropriate generalisation to the charged case is not clear 
to us at present. However we set up the full set of equations of motion, 
corresponding to a sufficiently general metric applicable to the charged case 
also. We then present the solutions for the uncharged case which reduce to 
the known ones in suitable limits and which may be thought of as decribing 
branes in a time dependent universe. They may also be thought of as de- 
scribing branes with a time dependent mass. We illustrate this possibility in 
a simplified context by making liberal assumptions. 

In string/M theory context, the static brane solutions are often inter- 
preted as coincident stacks of branes and antibranes [10]. In this context, see 
also [11] where more general static multiparameter solutions and their inter- 
pretation in terms of tachyon condensation are studied. The brane antibrane 
system has a tachyon mode indicating the instability due to brane antibrane 
pair annihilation. These pairs are expected to annihilate totally [12]. This 
time dependent decay process, also referred to as tachyon condensation, must 
have a gravity description atleast in the initial stages when the number of 
branes and antibranes is sufficiently large. It is then natural to hope that 
the present ansatz, suitably generalised, may describe such a process. We 
present a discussion of the necessary generalisations which, hopefully, are 
also sufficient. 

The organisation of this paper is as follows. In section 2 we present 
the equations of motion taking a sufficiently general metric. In section 3, 
we specialise to an ansatz similar to McVittie's, present the solutions and 
discuss various limits where they reduce to the known ones. In section 4 
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we discuss these solutions as describing branes with a time dependent mass. 
In section 5 we conclude with a discussion of the necessary generalisations 
of the present ansatz required to provide a gravity description of the brane 
decay process. 



2. Consider the following D— dimensional action for the metric Gmn, the 
scalar field (j), a (P + 2)— form gauge field strength F, and matter fields: 

S = — 2 j d D xV^G (p - - d M <P d M <p - 2(p - 2) , j + S mat (1) 

where G = det Gmn, A is a constant, 

I7i2 _ /~iMN rp2 r?2 rp rp M 2 -M P+2 

>' ( ' MN i MN — *MMi-M P+2 ' \ , 

and S mat is the action for matter fields which is assumed to be independent 
of 4> and F and has the energy momentum tensor Tmn given by 

5 = ~\j d D xV^G T MN 5G MN . 

The D— dimensional Newton's constant Gd is given by k 2 = 8ttGd- For 
D — 11, A = 0, and absent, the above action describes M theory branes, 
whereas for D = 10 and A = —1, +1, ^f- it describes fundamental strings, 
NS 5— branes, and DP— branes respectively. The equations of motion that 
follow from (1) can be written as 

Rmn = -ydM<t>dN(f> + S MN + k<2 Smn (2) 

V 2 = / e - F 2 (3) 
r 2(P + 2)! v ; 

d M (V^G e x * F MM *- Mp ^) = (4) 
V M T MN = (5) 
where the source terms S MN and Smn are given by 

S mn = 2(p + 2)\ + ^ ^ N ~~ P - 2 GmnF<2 ^) ( 6 ) 

Q _ rp GMN rp rp _ n M N rp ( r,\ 

• ' A I X 'MX Y) 2 ' ~~ +MN ■ \( ) 
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Consider the metric Gmn of the form given by 



ds 2 = 9ilv dx fl dx u + a (i) 5ij dy i dy j + S u ab {6 c ) d6 a d9 b (8) 

where fj, — 0, 1, • • • , (d — 1) ; i — 1, 2, • • • , m ; a = 1, 2, • • • , (n + 1) ; and g^ u , 
CT(j), and £ are functions of z' 1 only and are independent of y % and 9 a . Note 
that the dimensions of the spacetime is D = n + m + d+ I. The Ricci tensor 
Rmn for the above metric can be calculated straightforwardly and its non 
zero components can be written as 

Rij = ~ \ Gij V 2 In <7 (i) (9) 
R ab = R ab (uj) - i G ab V 2 In S (10) 

R,u = R,M - \ V 9 »Vl In (zS n+1 

_ 1 _ (w + 1) fl, , v 

4^U 2 Ao 4 ' 1 j 

In the above expressions, ( ) At = 3 /i ( ) , E = n™ i , is the covariant 
derivative with respect to g^ , R fiu (g) is the Ricci tensor of g^ , and R ao (uj) 
that of uj^ . Note that if u ab is the standard metric on a (n + 1)— dimensional 
sphere of unit radius, as will be the case here, then R a b(u) = n cu ab . Also, 
V 2 ?/> = G MN V mV 'n^P and, for any function ip(x^) which depends on x^ only, 
it is given by 

V 2 ^ = <T (v* V'tf + \ + ^±^j ^ . (12) 

For our purposes here, d — 2 and thus £) = n + m + 3, x M = (£, i?) 
where R is the isotropic radial coordinate, and y l are the coordinates of a 
m— dimensional compact space. Furthermore, we will take g^ u to be diagonal. 

Consider electric type P— branes with P < m. Then only F tR12 ...p (and 
related components obtained by permutations of indices) will be non zero 
and it depends on x^ only. Equation (4) then implies that 

FtR12 " P = ^_ g g gn+f ' 9 = det g^ = g tt g RR 
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where Q is a constant. It then follows that 



Q 2 e- 2X * 

i=P+l 



(P + 2)! " S ± 5"+ 1 ' S± n a d) 



and F^jv = GW^ 2 for MiV = (tt, RR,11, 22, PP) and = for 
other (MN) components. One then obtains that S MN = Gmn &(m) where 



V) " " 2 (D-2) v £l^J ' M -*» i2 » 1 » 2 »- 

^ = 2(^2) (f^) ' U ^ (13) 
Also, equation (3) becomes 



We will assume that the energy momentum tensor Tmn of the matter 
fields is of the perfect fluid form and is given by 

T t t = -p, T i j = 5 i jP(i) , T R R =p, T\ = 5\p. 

It then follows that Smn = Gmn S(m) where 

S(t) = ~d=2 (( D " 3 )' 9 + n ) ( 15 ) 
<%) = (p+(D-2)p (i) -n) (16) 

S(R) = S (a) = T ^ (p+ (D-2)p-U) (17) 



D-2 



with 



n = ^p w + (n + 2)p. 



3. In [3], McVittie had obtained solutions describing black holes in an 
FRW universe. See [4, 6, 5] for extensive discussions of these solutions, and 
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[6, 7, 8, 9] for some further generalisations. Here, we will study similar brane 
solutions in a time dependent universe. 

We consider here an ansatz similar to that of McVittie solutions and their 
higher dimensional generalisations [9]. The resulting solutions can be thought 
of as describing uncharged black branes in a time dependent universe with 
perfect fluid matter. 1 In this ansatz, F = , = a 7 , and 

ds 2 = -Zdt 2 + £ a 2 ^ (dy*) + a 2 ^\l + Y)» (dR 2 + R 2 dQ 2 n+1 ) (18) 
i=i 

where a(t) is a function of t only, 

Mlr^) 2 • (l9 » 

eti,(3, 7, and i? are constants, and rfQ n+ i is the standard line element on an 
unit (n + 1)— dimensional sphere. 

The above ansatz is chosen so that for a = 1 it describes static uncharged 
black brane solutions [1]. 2 Also, for R = the ansatz (18) describes an 
FRW universe where a(t) is the scale factor. For m = (3 = and R ^ 0, it 
describes McVittie solution in (n + 3) dimensions [9]. 

The line element in (18) can be written succinctly in terms of a function 
r(R, t) defined by 

r(R,t) = a R (1 +Y)$ . (20) 
It is then easy to verify that 

Z = 1 - -£ , r » = 4j R« (21) 

1 McVittie type solutions for Reissner-Nordstrom black holes in arbitrary dimensions 
are obtained in [7, 9]. Similarly, it should be possible to generalise the ansatz (18) and 
obtain solutions describing charged black branes in a time dependent universe with perfect 
fluid matter. The ansatz in (8) can indeed include the charged case but is far too general, 
and the resulting equations of motion are difficult to solve directly. Hence an ansatz in 
between (8) and (18) is needed. However, a few straightforward ones we tried did not 
work and the appropriate one is not clear to us at present. Hence we will consider only 
the uncharged case here. 

2 One can include arbitrary powers of Z in various terms in (18) and in which, for 
o=l, will then describe the multiparameter uncharged solutions given in [11]. 
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and that 



Ii = ftVI, ^ = ^ ; * =hXl f = £ (22) 

r r it z Z R 

where we have defined 

n _ d() _d() _a t 1 nr« 

The ansatz (18) now becomes 



m „ //7F? 2 \ 

ds 2 = -Zdt 2 + j:a 2a >(dy l ) )\ a 2 ? r 2 1^ + d{ll +1 j , = a? (23) 



where r(R,t) and Z (r(R,t)) are given in equations (20) and (21). 
One can now obtain the equations of motion. Defining 



m 

II 

i=l 



m 

= ^ «i , u = u + (n + 2) (/3 + v 7 ^) 



the equations of motion (2) and (3) can be written, after a long but straight- 
forward calculation, as 

hX(u+(n + l)p) = 
l(h t + (u-^)h 2 ^ = 

f(V(--f)//( = K 2 S { 



^(h t + (u-^)h 2 y^= (h t + uh 2 ) = k 2 S { 



(0 

(a) 



Y (E^ + (^ + 2)(^ + v / Z) 2 + ^-^ = -2« 2 p (24) 

where >% ia ) are given by equations (16), (17) and we have used F = and 
equations (9) — (11). 3 Equations (5) then follow from equations (24) by 
Bianchi identities. 



3 In (24), the second equation is obtained from equation (3) and the remaining ones 
from (2): the first, third, and fourth from MN = (tR, ij, ab or RR) in (2) and the last 
one is a combination of these and (2) with MN = tt. 
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Using equations (2) and (15) — (17), the Ricci scalar can be written as 

1 / m \ 

G MN R MN = - d M <j> d M 4> + k 2 \S it) + £ S {1) + (n + 2)S (a) j 

= \ d M cb d M <P + 2k 2 ^ (25) 

and can be expressed explicitly in terms of the quantities given in equation 

(23) by using equations (24) and noting that 

p - n m 

= E%) + ( n + 2 )^) -p- 

It follows from this explicit expression that if r 7^ then, generically, the 
Ricci scalar diverges as Z — > 0. For ctj = (3 = 7 = this divergence is 
absent if h t = 0, whereas for generic values of (aj,/3, 7) it is absent only if 
h t — h — 0. See [5] for a detailed study of this singularity for the D — 4, 
n — \ case. 

We now consider various limits where equations (24) reduce to the known 
ones. 

(I) For p = pi = p = and h — 0, equivalently a = constant, equations (24) 
are all satisfied. This is just the standard static black brane solution. 

(II) Let r = 0. Then Z = 1, X = and u = u + (n + 2) (/? + 1). Equations 

(24) become 



^(h t + uh 2 ) = 
aij (ht + u h 2 ^ = K 2 S(i) 
((3+1) (h t + uh 2 ) = n 2 S {a) 



h 2 (^a 2 + (n + 2)((3 + l) 2 + ^-u 2 ^ = -2k 2 p. (26) 

(II a) Standard FRW equations are obtained by setting 7 = 0, pi = p and, 
with no loss of generality, = 1 and (3 = in equations (26). Then 

u = (D-l) , U=(D-l)p, S {i) = S {a) =<^ 
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« (p - p) 



and we get 

+ P - 1) h D — 2 

(D -1){D -2)h 2 = 2k 2 P . (27) 

(II b) When p = pi = p = and a oat, Kasner type vacuum solutions can 
be obtained from equations (26) by setting n + 2 = formally and ignoring 
the 5(a)— equation. Then m = D — 1 and equations (26) imply that 

D-1 D-1 2 

«=E a i = 1 • E «* 2 + t = 1 • 
i=i i=i z 

The isotropic vacuum solutions studied in [13] are a special case of Kasner 
type solutions where ctj = a. Then 

1 9 2(D-2) 

(III) Let r 7^ and h ^ 0. Equations (24) then imply that, generically, 

u + (n + 1)13 = 7 = . 

If = also, equivalently a, = 0, for i = 1,2, ■ ■ ■ ,m, then it = and, 
hence, (3 = 0. Then, as follow from equations (24), 

kt + (n + 2)/, 2 = ^S ia) = K2{p - p) 



v ' w n+1 

(n + l)(n + 2)/i 2 = 2/€ 2 p 

which imply that p is homogeneous and is a function of t only. Note that 
if m — 0, ie. if there are no compact directions, then in equations (24) 
a (i) — Pi — formally and <%)— equations are to be ignored. Then (3 — 
and equations (24) reduce to the above ones with D = n + 3. These are just 
the higher dimensional, spatially flat, McVittie solutions given in [9]. 

4. When o = 1 the function Z in (21) is time independent and the ansatz 
(18), equivalently (23), describes uncharged branes with ADM mass M oc r$ 
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[11]. When a(t) is time dependent, the function Z is also time dependent 
and then the ansatz (18) can naturally be thought of as describing uncharged 
branes with time dependent mass, the time dependence being dictated by 
a(t). 

In string/M theory context the brane solutions (18) with a = 1 are often 
interpreted as coincident stacks of equal number, N, of branes and antibranes 
with N ^> 1 [10]. 4 Such a system has a tachyon mode indicating the 
instability due to brane antibrane pair annihilation. As these pairs decay, 
the tachyon field rolls down from the maximum of its potential and condenses 
at the potential minimum. At the end of this process the brane antibrane 
pairs would have annihilated totally and disappeared. 5 

In recent years, enormous progress has been made in understanding the 
tachyon condensation process in the case of a single brane antibrane pair. 
See the recent comprehensive review [12] and references therein. However, 
not much is known in the case of stacks of brane antibrane pairs with N ^> 1. 
Nevertheless, it is reasonable to assume that the dynamics is qualitatively 
similar to that of a single brane antibrane pair; namely, that all the brane 
antibrane pairs will annihilate totally at the end. 

For N sufficiently large, the above dynamic process must have a gravity 
description. During this process N, and the mass M , will decrease in time. 
Such a description must be valid atleast in the initial stages when N(t) 
remains sufficiently large. Then it is natural to hope that an ansatz similar 
to (18) may describe such a process. 

That the ansatz (18) can describe time dependent mass when a(t) is time 
dependent can be shown by using the quasilocal Hawking mass [5, 14]. Here, 
for the purpose of illustration, let us instead imagine that the system remains 
static for t < with a(0) = 1, that the time dependence is switched on at 
time t = 0+, and off at t = t _ > with a(t ) = a , and that the system 
remains static again for t > t . 6 Let us also assume that a ^ 1 and that the 

4 More general static multiparameter solutions and their interpretation in terms of 
tachyon condensation are studied in [11]. 

5 In the charged case the numbers ./V and N of branes and antibranes are not equal and 
so the decay process will end in extremal branes with charge oc (N — N). 

6 In the case of tachyon condensation, switching on the system is analogous to nudging 
the tachyon field away from its potential maximum towards the minimum; switching off, 
while its meaning is not entirely clear, is perhaps analogous to the condensation at the 
potential minimum if it happens in a finite time. 
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parameters in the ansatz (18) obey u + (n + 1)(3 = for t > to, assumptions 
which seem physically reasonable. 

We then have (a,h,h t ) = (1,0,0) for t < and = (a ,0,0) for t > t . 
The initial and final masses, denoted as M and M respectively, can then 
be calculated 7 in the standard way [11] without requiring the concept of 
quasilocal mass [5, 14]. They are given by 

M= ( " +1) ^+ lKnr ° , M = aT n "M=^ (28) 

where uj n+ i is the area of an unit (n+ 1)— dimensional sphere, V m is the initial 
volume of the m— dimensional compact space, and the relation u + (n+ 1)(3 = 
is used. So, the final mass M is different from the initial mass Mif oj^ 1. 
Thus if Oq > 1 then M < M indicating that the final mass decreases. Note 
that if u = 0, which is necessarily the case if m = 0, then (3 = and M = M. 
This is the case for McVittie solutions since m = u = there [3, 8, 9]. 

Bekenstein-Hawking entropies Sbh an d Sobh, for t < and t > to re- 
spectively, can also be calculated and are given by 



2vr cu n+1 V m r " +1 

~>BH — 



K 2 



Sobh — ao +<n+1 ' )/3 Sbh — Sbh ■ (29) 



Thus the final entropy S bh and initial entropy Sbh are equal here which is a 
consequence of the relation u+ (n + 1) (3 = which, in turn, is a consequence 
of the absence of the source term S t R as can be seen from (24). This indicates 
that a non zero source term StR is likely to be necessary to model a system 
such as decaying brane antibrane pairs where entropy decrease is expected 
on physical grounds. However, a non zero S t R alone may not be sufficient, 
see the discussion below in section 5. 



5. Clearly a more general ansatz is needed to model brane antibrane decay 
process, equivalently tachyon condensation. We now discuss the necessary 
generalisations which, hopefully, are also sufficient. 

In the decay of a single brane antibrane pair (N = 1), the tachyon field 
resides in the brane world volume. Also, in the weak coupling limit where 
the string coupling constant g s — > 0, the total energy on the brane remains 

7 The calculation of M is straightforward once one defines Rfi n — a\ +l3 R in (18). 
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constant as the tachyon rolls down from the potential maximum to the mini- 
mum. This is because the emission of closed string radiation is suppressed in 
the limit g s — > 0. However, for non zero g s , it is expected that closed strings 
are radiated into the transverse non compact space, thereby decreasing the 
total energy on the brane. For more details see [12] and references therein. 

Assuming that similar qualitative features persist also for the case when 
N ^> 1 and g s is non vanishing then suggests the following generalisations: 

• A delta function source S\\ 5{R trans ) must be included with S\\ given by 
the energy momentum tensor of a world volume scalar field T rolling down 
its potential V(T). 8 Note that, even for the N = 1, g s — > case, V(T) is 
not known from first principles. For the N ^> I and non vanishing g s case 
required here, one has to first find V(T) from first principles or, as in the 
N — 1 case, find an empirical form with correct properties. 

Instead of a delta function source one may include, more generally, a 
smeared source localised close to the brane world volume and assume different 
metrics, e.g. the metric in (8) with (g^, cr (i) , S) in and (g^, cr (i) , S) out , and 
different source terms in this region and in the region outside it. Appropriate 
junction conditions are then to be imposed to ensure continuity across the 
boundary. See [4] for an example. 

• Source terms representing outgoing radiation in the transverse space must 
be included. Such terms will have non zero S t R. It may then be possible to 
have Sqbh < Sbh, see equation (29) and related comments. 

• Area increase theorems of General Relativity [15] forbid the decrease of 
horizon area and hence of Bekenstein-Hawking entropy Sbh- Now, if the 
brane decay dynamics is to result in the decrease of horizon area, and thus of 
Sbh, as seems physicaly reasonable then these theorems must be invalidated, 
presumably by violating some energy conditions assumed in their proofs. 

It is not clear if the outgoing radiation terms, mentioned above, will 
suffice. So, in general, one may include empirically a perfect fluid source 
with an appropriate equation of state so that some combination of outgoing 
radiation and this perfect fluid may invalidate the area increase theorems and 
result in the decrease of Sbh- In this context note that black hole mass and, 
consequently, entropy is known to decrease due to the accretion of phantom 

8 The source term and the energy momentum tensor are related as in equation (7). 
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energy (a perfect fluid with - < — 1) which violates energy conditions, see [16] 
and references therein. It will be quite interesting if similar energy source (s) 
were to originate in string theory in the context of brane decay. 

In conclusion, we note that incorporating the above generalisations in the 
ansatz for the metric and the source terms Smn is comparatively easy. Indeed 
the ansatz given in equation (8) for the metric is quite general and suffices 
also when the above generalisations are incorporated. However, it is a tall 
order to solve the resulting equations of motion, extract the salient features, 
and interpret them in terms of the dynamics of brane decay, equivalently 
tachyon condensation, in string theory. Nevertheless such a project seems 
worthwhile in view of the insights it is likely to offer. 
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